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The PHENIX experiment at the Relativistic Heavy Ion Collider has measured (p meson production 
and its nuclear modification in asymmetric Cu+Au heavy-ion collisions at ,/snn = 200 GeV at both 
forward Cu-going direction (1.2 < y < 2.2) and backward Au-going direction (—2.2 < y < —1.2), 
rapidities. The measurements are performed via the dimuon decay channel and reported as a 
function of the number of participating nucleons, rapidity, and transverse momentum. In the most 
central events, 0%-20% centrality, the (p meson yield integrated over 1 < pt < 5 GeV/c prefers a 
smaller value, which means a larger nuclear modification, in the Cu-going direction compared to 
the Au-going direction. Additionally, the nuclear-modification factor in Cu-|-Au collisions averaged 
over all centrality is measured to be similar to the previous PHENIX result in d-|-Au collisions for 
these rapidities. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

The Relativistic Heavy Ion Collider (RHIC) acceler¬ 
ator and its four experiments have previously provided 
extensive experimental evidence to confirm the forma¬ 
tion of a deconfined state of nuclear matter, referred to 
as the quark-gluon plasma (QGIH, in the initial stages of 
high-energy heavy-ion collisions [ll-lj . Currently, a major 


* Deceased 
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objective in the field of high-energy nuclear physics is to 
characterize the properties of the QGP in a quantitative 
way. The (p meson is a useful probe for studying the QGP 
properties, because it is sensitive to several aspects of the 
collision, including modifications of strangeness produc¬ 
tion in bulk matter Due to its small inelastic cross 

section for interaction with nonstrange hadrons [^, |^, 
the (p meson is less affected by late hadronic rescattering 
and may reflect the initial evolution of the system. Be¬ 
ing composed of a nearly pure strange anti-strange (ss) 
state, the cp meson puts additional constraints on models 
of quark recombination in the QGP. 

The study of the QGP typically involves compar- 
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isons of different observables measured in nucleus-nucleus 
(A-l-B) collisions and in proton-proton (p+p) collisions 
at the same center-of-mass energy. Modifications in the 
A-l-B collisions with respect to p+p collisions could be in¬ 
terpreted as being due to the hot nuclear matter (HNM) 
- possibly QGP - being produced. However, nuclear 
modifications could be present in the initial state of the 
collisions even if no QGP is produced. These effects, 
typically referred to as cold nuclear matter (CNM), may 
include the modification of parton distribution functions 
(PDFs) in a nucleus Q, initial-state energy loss [13 , and 
the Cronin effect, which is often attributed to multiple 
scatte ring of the incoming parton inside the target nu¬ 
cleus 11, ill- CNM effects can be probed with d-|-Au 
collisions. PHENIX has previously measured 4> meson 
production in d-l-Au collisions at forward, mid- and back¬ 
ward rapidities [l3| . Suppression was observed in the for¬ 
ward (d-going) direction, where small-x partons from the 
Au nucleus are probed, and an enhancement was seen in 
the backward (Au-going) direction. Similar behavior was 
previously observed for inclusive charged hadrons and 
open heavy flavor in d-|-Au collisions 1J,[I5|, potentially 
indicating similar particle production and modification 
mechanisms. 


The rapidity dependence y of particle production in 
asymmetric collisions with a smaller-A projectile and 
a large-A target, provides a way to investigate both 
hot and cold nuclear-matter effects. Previous J/4’ me¬ 
son data in Cu-fAu collisions M showed that the ra¬ 
tio of forward (1.2 < y < 2.2, or Cu-going) to back¬ 
ward (—2.2 < y < —1.2, or Au-going) J/'0 modification 
was comparable in both sign and magnitude to that ex¬ 
pected from GNM effects. The (f) meson is composed 
of lighter closed flavor (ss) and its production from 1.0 
GeV/c to 5.0 GeV/c involves a mix of soft and hard pro¬ 
cesses and would provide a link between heavy flavor and 
lighter mesons. Gomparison of the (j) meson production 
in Gu-|-Au and d-|-Au systems and to J/ip production in 
Gu-|-Au collisions may shed light on the mixture of HNM 
and CNM effects on ())-meson production. 

The production of (p mesons has already been mea¬ 
sured at PHENIX in P +P. d-fAu, Cu-fCu, and Au-fAu 
at midrapidity |17Hl9l | and in p+p and d-|-Au at forward 
and backward rapidities [H, l20|| over a wide range in 
Pt- Previous measurements from Au-fAu and Cu-fCu 
collisions in a similar momentum range were found 
to be consistent with HNM effects and exhibited large 
flow anisotropies. The STAR Collaboration has also pre¬ 
viously measured (p meson production at midrapidity in 
Cu-I-Cu and Au-fAu collisions [2l|,[2^. p meson produc¬ 
tion has also been measured by the ALICE Collaboration 
at large rapidity in p+p and p-bPb collisions [2^ and at 
midrapidity in Pb-|-Pb collisions (^ . 

In this paper, the production of p mesons is determined 
at forward and backward rapidities via dimuons recon¬ 
structed in the PHENIX muon spectrometers in Cu-|-Au 
collisions at ^/spppp= 200 GeV recorded in 2012. The par¬ 
ticle multiplicity at these rapidities in heavy-ion collisions 


results in large combinatorial backgrounds and produces 
a challenging environment for p meson measurements. 
Previous measurements were thus limited to smaller col¬ 
lision species. A procedure for removing the background 
is detailed and a measurement of the p meson nuclear 
modification factor i?cuAu in Cu-f Au collisions at for¬ 
ward and backward rapidities is presented versus y, pr, 
and the number of participating nucleons. 


II. EXPERIMENTAL SETUP 


The PHENIX detector is described in detail in [25l| . 
and a schematic of the 2012 setup is shown in Fig. [T] 
This analysis uses the dimuon decay channel of the p 
meson. The detectors relevant for this measurement 
are forward and backward muon spectrometers [2^, the 
two beam-beam counters (BBCs) [^, the silicon vertex 


tracker fVT Xl I2g 
tor (FVTX) m 


and the forward silicon vertex detec- 



18.5 m= 60 n 


Side View North 


L. 

ZDC North 
MuID 


FIG. 1. (color online) The 2012 setup of the PHENIX detec¬ 
tor. 


This study used minimum bias (MB) events triggered 
by the BBGs. The BBCs comprise two arrays of 64 
Cerenkov counters covering the pseudorapidity range 
3.1 < \r]\ < 3.9. The MB trigger required two or 
more counters firing on each side and a z-vertex selec¬ 
tion around the nominal center of the detector accep¬ 
tance [l^. The MB trigger fired on 93±3% of the 5.2±0.2 
b total inelastic Cu-|-Au cross section. In this case, the 
z-vertex was measured by the BBCs with a resolution of 
cr 2 ~ 0 .5-2.0 cm, depending on the event multiplicity. 

The collision point is determined in x, y and z by the 
two vertex detectors, VTX and FVTX, with a resolution 
of better than 100 microns. The VTX and FVTX detec¬ 
tors were installed in 2011 and 2012 to provide precise 
particle vertexing and tracking in the central and for¬ 
ward/backward rapidities. Covering approximately the 
same rapidity range as the existing muon spectrometers, 
the FVTX is composed of two endcaps, each with four 
stations that are perpendicular to the beamline and com¬ 
posed of silicon mini-strip sensors that have a 75 micron 
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pitch in the radial direction and lengths in the (j) direction 
varying from 3.4 mm to 11.5 mm. The VTX, which sur¬ 
rounds the collision region at PHENIX, comprises four 
layers of silicon sensors. The inner two layers and outer 
two layers are composed of 30 pixel ladders and 44 strip- 
ixel ladders, respectively. 

The muon system is separated into the north and south 
muon arms. Each arm comprises four subcomponents: 
an absorber material, a magnet, a muon tracker (MuTr), 
and a muon identifier (MulD). Initially, the absorbers 
were composed of 19 cm copper and 60 cm iron, but an 
additional 36.2 cm of stainless steel was added in 2010 
to help decrease the hadronic background. Eollowing the 
absorber in each muon arm is the MuTr, which comprises 
three sets of cathode strip chambers in a radial magnetic 
field with an integrated bending power of 0.8 T-m. The 
final component is the MulD, which comprises five alter¬ 
nating steel absorbers and larocci tubes to further reduce 
the number of punch-through hadrons that can be mis¬ 
takenly identified as muons. The backplates of the mag¬ 
nets provide the first absorber layer for the muon identi¬ 
fier systems. The backplate of the south muon magnet is 
10 cm shorter than the backplate of the north muon mag¬ 
net, resulting in less total absorber material in the south 
arm than the north arm, and thus a slightly different 
momentum acceptance. The muon spectrometers cover 
the pseudorapidity range 1.2 < \r]\ < 2.2 over the full az¬ 
imuth. Muon candidates are identified by reconstructed 
tracks in the MuTr matched to MulD tracks, where at 
least one of the tracks from a pair of muon candidates 
in the same event penetrates through to the last MulD 
plane. The minimum momentum needed for a muon to 
reach the last MulD plane is ^3 GeV/c. 


III. DATA ANALYSIS 
A. Dataset and quality cuts 

In this analysis, (j) meson candidates are selected from 
two reconstructed muons in the RHIC Cu-|-Au dataset 
from 2012. The (j) meson invariant yields are then mea¬ 
sured and used to calculate the nuclear modification fac¬ 
tor i?cuAu, which is compared to results from other sys¬ 
tems. For this analysis, 4.73 billion (£ = 0.97 nb“^) 
sampled MB events were used within ±10 cm z-vertex 
and 0%-93% centrality. The total inelastic cross section 
for Cu±Au collisions at 200 GeV was estimated by a 
Glauber simulation to be 5.2±0.2 b. 

A set of quality assurance cuts is applied to the data 
to select good muon candidates and improve the signal- 
to-background ratio. These cuts are summarized in Ta¬ 
ble n The collision z-vertex is required to be within 
±10 cm of the center of the interaction region along the 
beam direction, as measured with the BBCs. The MuTr 
tracks are required to match the MulD tracks at the first 
MulD layer in both position and angle. In addition, only 
dimuon candidates in which at least one track penetrated 


to the final MulD layer are selected. Furthermore, the 
track is required to have greater than a minimum number 
of possible hits in the MuTr and MulD, and a maximum 
allowed is applied to both the track and vertex deter¬ 
mination. There is a minimum allowed single muon mo¬ 
mentum along the beam axis, Pz, which is reconstructed 
and energy-loss corrected at the collision vertex. Finally, 
this analysis is restricted to the dimuon px range of 1 — 5 
GeV/c. This limitation is due to the large backgrounds 
and small acceptance at low px and small statistics at 
high px, preventing signal extraction of the (f> meson. 
The events are sorted into centrality classes using the 
combined charge from both BBCs [2. The number of 
binary collisions Vcoii and number of participating nucle¬ 
ons Apart are extracted from a Glauber simulation [l^ . 

B. Background subtraction 

The PHENIX muon spectrometers have a small ac¬ 
ceptance for (j) mesons. Going from the most peripheral 
centrality bin, 40%-93%, to the most central bin, 0%- 
20%, the signal-to-background ratio decreases from 0.28 
to 0.067 in the Cu-going direction (1.2 < y < 2.2) and 
from 0.37 to 0.090 in the Au-going direction (—2.2 < y < 
— 1.2). Due to the very low signal-to-background ratio, 
particularly in the most central events, the background 
subtraction is of crucial importance. Accordingly, several 
different background subtraction methods were explored 
and compared. 

The invariant mass distribution is formed by com¬ 
bining muon candidate tracks of opposite charge. This 
unlike-sign invariant mass spectrum contains the (j), p 
and w mesons as well as both uncorrelated and corre¬ 
lated backgrounds. The uncorrelated backgrounds come 
from random combinatorial associations of muon candi¬ 
dates, while the correlated backgrounds arise from open 
charm decay (e.g., DD where both decay semileptoni- 
cally to muons), open beauty decay, rj meson and w me¬ 
son Dalitz decays and the Drell-Yan process. These cor¬ 
related backgrounds are described in Sec. IIIC. The un¬ 
correlated combinatorial background is accounted for via 
two methods: (1) like-sign dimuons and (2) event mixing. 

First, the uncorrelated combinatorial background is 
estimated through the like-sign background subtraction 
technique, which is generally associated with the assump¬ 
tion that the like-sign dimuon pairs come purely from 
combinatorial processes without any correlation between 
muons. It follows that the like-sign distribution can be 
subtracted from the unlike-sign distribution according to 
the relationship described in Eq. [1] 

A+_=FG+_-FG±±, (1) 

where A_|_is the uncorrelated background subtracted 

signal and FG^ _and FG±± are the unlike-sign and like- 

sign dimuon pairs, respectively, corresponding to pairs 
formed within the same event. The like-sign distribution 
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TABLE 1. Quality cuts for 0 meson signal extraction in Cu+Au collisions. 


Variable 

Au-going 

Gu-going 

Meaning 

Zvtx (cm) 

< 10 

< 10 

Collision vertex along the beam direction as measured by the BBCs 

pDGO 

< 90 

< 50 

Track momentum times the spatial difference between 

(GeV/c- cm) 



the MuTr track and MuID track at the first MuID layer 

pDDGO 
(GeV/c- radian) 

< 30 

< 45 

Track momentum times the slope difference between 
the MuTr track and MuID track at the first MuID layer 

Track 

< 5 

< 10 

X^/NDF of the p track 

Lastgap 

one track > 2 

other track > 4 

one track > 2 

other track > 4 

Last MuID plane that the p track penetrated 

nidhits 

ntrhits 

Xvtx 

Dimuon pr (GeV/c) 
\pz\{GeV/c) 

>(2xlastgap —1) >(2xlastgap —1) Number of hits in the MuID, out of the maximum 10 
>11 >10 Number of hits in the MuTr, out of the maximum 16 

<4 <7 x^/NDF of the dimuon track with the vertex 

1 — 5 1 — 5 Transverse momentum of the dimuon pair 

> 2.4 > 2.5 Momentum of the fi along the beam axis 


FG±± is normalized to a quantity that is more precise 
and not sensitive to differences in the detector acceptance 
between like-sign and unlike-sign pairs. This background 
normalization is described in Eq. 


2Jj FG++dm J FG—dm 

FG,, = iFG,, + re__) ’ 

( 2 ) 

where m is the dimuon invariant mass, and the integra¬ 
tion is carried out in the range 0.2 < m < 5.0 GeV/c^. 

In parallel to the like-sign technique, the uncorrelated 
background is also estimated through the event mixing 
technique. In the standard event mixing method, muons 
from different events are randomly associated to produce 
a background distribution of uncorrelated dimuon pairs. 
Events were mixed with partners from within the same 
2%-centrality and 1-cm z-vertex bins in order to minimize 
the systematic uncertainties. The mixed-event back¬ 
ground distributions {BG) were generated with about 8 
times higher statistics than the actual background and 
then normalized to match the same-event foreground 
{FG). The normalization factor also accounts for slightly 
different multiplicities from mixing of slightly different 
events. Although a mass-dependent technique was devel¬ 
oped for this analysis, a standard event mixing technique 
is described in advance. In previous PHENIX analyses, 
the normalization factor a was calculated as described in 
Eq.[3] 


/ FG++dm / FG -.dm 
f BG-i-+dm f BG _dm’ 

where E'G_|_+ and FG _are the like-sign pairs from the 

same event and i?G_|_+ and BG _are the like-sign pairs 

from mixed events. 

After subtracting and fitting the resonances as well 
as the remaining correlated background, the yields from 
mixed-event background subtraction are consistent with 
the yields from the like-sign technique within statisti¬ 
cal uncertainties. The event mixing technique is used 
in this analysis due to the statistical limitations of the 
like-sign technique. The differences between the like-sign 
and event mixing techniques are used to determine one 
component of the systematic uncertainty on the yield, as 
described later in Sec. IIIF. 

In this method, each term in the square root of 
Eq.l was integrated over all mass, introducing a mass- 
independent normalization factor [16l . |3l| . Dimuons from 
same events are less likely to be reconstructed in close 
proximity to each other than those in mixed events, 
resulting in a larger relative number of mixed-event 
dimuons at low mass, where the opening angle is small, 
than at higher mass. Therefore, the normalization fac¬ 
tor, which is simply a ratio of the like-sign same-event 
dimuons to like-sign mixed-event dimuons, drops at lower 
masses. Because this normalization factor depends on 
mass, particularly in the cj) meson region, it became nec¬ 
essary to introduce a mass-dependent normalization, as 
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FIG. 2. (color online) The event mixing normalization factor a versus mass. This factor shows a dependence on mass, 
particularly in the low mass region. The error function, which was used to fit a, can also be seen in the plot along with the fit 
parameters and goodness of fit. 


described in Eq. |4l rather than the more commonly used 
mass-integrated normalization from Eq. [3l 


a(m) 


! FG++{m)FG—{m) 
BG++ (m)BG _(m) 


(4) 


This mass-dependent normalization factor is then fit 
as a function of mass, and the fit function - rather than 
the integrated normalization factor - is multiplied to the 
unlike-sign mixed-event background to get the normal¬ 
ized background spectrum yjG'normahzed^ 


^^„ormalized(^) = a{m) X BG+-{m). (5) 

Several fitting functions were tested, including a poly¬ 
nomial and an error function. The error function, which 
is used in the final analysis, is described in Eq. (SJ where 
girn) is the error function and po, pi and P 2 are free pa¬ 
rameters of the fit. A plot of the normalization factor as 
a function of mass fit with an error function is shown in 
Fig.H 

g{m) =poX Fvii———) (6) 

P2 

The application of event mixing to describe and sub¬ 
tract backgrounds in the 4> meson mass region is shown 
in Fig. [3l where the open squares represent the mixed- 
event background and the closed circles are the unlike- 
sign spectrum. Before background subtraction, the p+uj, 
(j) and J/ip peaks are clearly seen. 


C. Signal extraction and correlated background 

After the mixed-event background subtraction, there 
is still some correlated background remaining. In previ¬ 
ous PHENIX analyses, it was shown that heavy flavor 
(charm and beauty) contributions were negligible in the 
(j) meson mass region for p+p and d-pAu collisions at 
200 GeV [H, [13. Simulation studies showed that rj me¬ 
son Dalitz decays are one possible contributor to the cor¬ 
related background. The correlated background is well 
described by the function in Eq. [7] 


/(m) = exp(a • m) + b + c ■ m, (7) 

where a, b and c are free parameters of the fit f{m). 
Accordingly, the correlated background in real data are 
also fit with the function described in Eq. [71 as shown 
in Fig. 21 where the mass distribution after mixed-event 
background subtraction is shown. Several other fit func¬ 
tions and fit ranges were tested and used to estimate a 
systematic uncertainty. 

The (j) and uj meson signals are each described by a 
Gaussian and the signal from the p meson by a Breit- 
Wigner distribution, as shown in Fig. 21 along with the 
correlated background description. The (j) meson mass 
resolution is ^90 MeV/c^. The PHENIX muon arms are 
not able to resolve the p and ui peaks separately, so a 
combined fit is made. All fit parameters are constrained 
but allowed to vary, except the ratio of the yield of p 
mesons to that of p + uj, which is set as a constant based 
on the expected ratio between their cross sections and 
branching ratios. The data are binned as a function of 
Pt, y and centrality over the range 1 < pt < 5 GeV/c, 
1.2 < \y\ < 2.2, and 0%-93% centrality. 
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FIG. 3. (color online) The unlike-sign spectra and combinatorial background described with event mixing for 1.2 < y < 2.2 
(Cn-going direction) and —2.2 < y < —1.2 (An-going direction). The p + ui, (j) and J/ip peaks are clearly visible before 
background subtraction. The mass bin width is 71 MeV as marked on the vertical axis. 



(GeV/c^) 


FIG. 4. (color online) The dimuon mass spectra for 1.2 < y 
direction) after subtracting mixed events and fitting the ij) and 
mass bin width is 71 MeV as marked on the vertical axis. 


D. Detector acceptance and reconstruction 
efficiency 

The product of detector acceptance and reconstruction 
efficiency, Asrec, of dimuon decays of (j) mesons is deter¬ 
mined by the full event reconstruction of the (j) meson 
signal obtained from pythia 6.42 [s^, run through a full 
GEANT3 [13 simulation of the 2012 PHENIX detector 
setup, and embedded in the MB real-data background. 
The embedded simulated events are then reconstructed 
in the same manner as data with the same cuts applied 



C 2.2 (Gu-going direction) and —2.2 < y < —1.2 (Au-going 
p + oj peaks and the remaining correlated background. The 


as in the real data analysis. The background subtraction 
and signal extraction are also handled in the exact same 
manner as in real data. The Ae^^c is then calculated as 
the number of reconstructed (j) meson candidates divided 
by the number of (p mesons generated in pythia, both 
within an appropriate kinematic bin. As previously men¬ 
tioned, the south arm has a smaller amount of absorber 
material, causing a larger acceptance in the south arm 
(Au-going direction) than in the north arm (Cu-going 
direction). In addition, the Asrec bas a centrality and 
Pt dependence. Specifically, for the lower pr bin (1-2.5 
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GeV/c), Affrec = 1-21 X 10“^ in the Cu-going direction 
and 1.86 x 10“^ in the Au-going direction, while for the 
higher pT bin (2.5-5 GeV/c), As^ec = 1-69 x 10“^ in 
the Cu-going direction and 1.81 x 10“^ in the Au-going 
direction. The centrality dependence is not as strong, 
with the values going from Asrec = 2.23 x 10“^ in the 
Cu-going direction and 2.37 x 10“^ in the Au-going di¬ 
rection at 0%-20% centrality to Ae-cec = 2.41 x 10“^ in 
the Cu-going direction and 3.83 x 10“^ in the Au-going 
direction at 40%-93% centrality. 


E. Invariant yields and nuclear modification factors 

The invariant yield is calculated according to the rela¬ 
tion: 


BR 


dydpT 


1 N 

^y^PT As YecAfevt 


( 8 ) 


where BR is the branching ratio to dimuons {BR((j) —>■ 
pApi~) = (2.89 ± 0.19) X 10“^ H): -^evt is the number 
of sampled MB events within the relevant centrality se¬ 
lection (Aevt = 4.73 X 10® for the 0%-93% selection), N 
is the number of observed (j) mesons, and Ay and Apx 
are the bin widths in y and pt, respectively. To eval¬ 
uate the nuclear matter effects on (j) meson production 
in Cu-|-Au collisions, the (f) meson yields in Cu-|-Au colli¬ 
sions are compared to those measured in p+p collisions at 
the same energy after scaling by the number of nucleon- 
nucleon collisions in the Cu-|-Au system, A^coii- This ra¬ 
tio is called the nuclear modification factor i?cuAu, and 
is defined as: 


T^CuAu — 


A^cuAu 

dydpT 


Acoll X 


d^Afpp ■ 

dydpT 


(9) 


The p+p reference data used in the RcuAu are from 
Ref. [20|. Because the rapidity and pr binning in the 
Cu-|-Au analysis differs from that in the p+p analysis, 
the p+p invariant yields were re-measured using the same 
binniM as the Cu-|-Au yields and in a manner similar to 
Ref. . The sampled luminosity of the p+p data used 
in this analysis corresponds to £ = 14.1 pb“^ [^ . 


F. Systematic uncertainties 

The systematic uncertainties associated with this mea¬ 
surement are categorized as Type-A, Type-B or Type-C. 
Type-A refers to point-to-point uncorrelated uncertain¬ 
ties that allow the data points to move independently 
with respect to one another. They are added in quadra¬ 
ture with the statistical uncertainties and represented on 
the plots as an error bar. Type-B uncertainties are cor¬ 
related point-to-point, which means the points move co¬ 
herently. All sources of Type-B uncertainty are added 


TABLE II. Systematic uncertainties included in 
yield calculations. 

the invariant 

Type 

Origin 

Value 

A 

Signal extraction 

2-31% 

B 

MuID efficiency 

2% 

B 

MuTr efficiency 

2% 

B 

■^^Tec 

13% 

B 

(j) candidate selection 

3% 

B 

Like-sign background subtraction 

5% 

C 

MB trigger 

3% 


TABLE III. Systematic uncertainties included in 
modification factor calculations. 

the nuclear 

Type 

Origin 

Value 

A 

Signal extraction 

2-31% 

A 

p+p reference (integrated centrality only) 

5-13% 

B 

MuID efficiency 

4% 

B 

MuTr efficiency 

2% 

B 

■^^rec 

13% 

B 

(p candidate selection 

3% 

B 

Like-sign background subtraction 

5% 

B 

Neon (centrality bins only) 

5-10% 

C 

MB trigger 

10% 

C 

Acoii (integrated centrality only) 

5% 

C 

p+p reference (centrality bins only) 

11% 


in quadrature and displayed as boxes around the data 
points. Finally, Type-C refers to the global uncertainties 
which allow the data points to move together by an iden¬ 
tical multiplicative factor. The Type-C uncertainties are 
given in the legends of the plots. 

Several systematic uncertainties are evaluated for this 
analysis. For the signal extraction uncertainty, different 
fits and parameters are tested for the background nor¬ 
malization factor, the correlated background, the p + lu 
signal, and the (j> meson signal. This is done separately 
for each kinematic bin, and a 2-31% systematic uncer¬ 
tainty is assigned, with the largest uncertainty on yields 
extracted from the most central events. This is because 
the high multiplicity in central collisions results in large 
combinatorial backgrounds and a very small signal-to- 
background ratio. It is important to note here that the 
signal extraction uncertainty was primarily dominated by 
the fluctuations in the correlated background. The p+p 
reference uncertainty comes from the uncertainty on the 
(j) yields in the p+p reference 0. There is a 4% system¬ 
atic uncertainty from the MuID efficiency and a 2% un- 
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certainty from the MuTr efficiency in p+p collisions . 
In Cu+Au collisions, the MuTr efficiency uncertainty re¬ 
mains the same, while the MuID efficiency uncertainty 
drops down to 2% [I^. For the Ae^ec uncertainty, the 
Pt and y distributions in pythia are changed to match 
the slope of the distributions in real data, and allowed 
to vary over the range of the error bars in data, yielding 
a 13% systematic uncertainty. Real data and simulation 
inconsistencies in each of the muon identification cuts 
listed in Table Hare also evaluated. They can affect the 
yields by 3%, which is assigned as a systematic uncer¬ 
tainty on the (j) meson candidate selection. The like-sign 
background subtraction uncertainty of 5% comes from 
differences in the yields when using the like-sign method 
or the event mixing method. The A^coii uncertainty of 
5-10% arises from the fact that Acoii carries a statistical 
uncertainty itself. Finally, the MB trigger efficiency un¬ 
certainty was 10% in the p+p reference [2^ and 3% in 
Cu-|-Au collisions [l^ . All of these systematic uncertain¬ 
ties are tabulated in Tables m and m 


IV. RESULTS 

The invariant yields for 1 < pT < 5 GeV /c (p mesons 
are calculated as a function of centrality, y and pr as 
described in Eq. [51 The results are summarized in Ta¬ 
bles |IV]-|Vll Similarly, the nuclear modification factors 
are formed from the invariant yields using Eq. [9] and tab¬ 
ulated in Tables IVlIl - IIXI 

Fig- El shows the invariant yield as a function of the 
number of participating nucleons Apart- In Fig. |6l the 
dependence of the invariant yield on transverse momen¬ 
tum Pt is shown. The invariant yield as a function of 
rapidity is plotted in Eig. 0 More p mesons are pro¬ 
duced in the Au-going direction (—2.2 < y < —1.2) than 
in the Cu-going direction(1.2 < y < 2.2). This may be 
explained by the larger multiplicity in the Au-going di¬ 
rection coupled with a mixture of both HNM and CNM 
effects. 

Although the invariant yields are interesting on their 
own, the nuclear modification factor is studied in order 
to evaluate the effects of hot and cold nuclear matter 
on (f> meson production in Cu-|-Au collisions at y/s^= 
200 GeV. 

The nuclear modification factor as a function of Apart 
is shown in Fig. |8l There is a dependence of i?cuAu 
on both centrality and rapidity. In the Au-going direc¬ 
tion, the i?cuAu is greater than unity for all centralities. 
The rapidity dependence is similar to the trend observed 
by PHENIX for p —>■ in d-|-Au collisions [l^ as 

well as measurements made by the ALICE Collabora¬ 
tion at large rapidity in p-|-Pb collisions at 5.02 TeV at 
the Large Hadron Collider [2^, where an enhancement 
was observed in the Pb-going direction while the p-going 
direction was either suppressed or consistent with unity 
depending on the pr range. 

To further understand the relative roles of different 
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FIG. 5. (color online) Invariant yield as a function of the 
number of participating nucleons for 1.2 < [i/l < 2.2 and 
1 < Pt < 5 GeV/c. The centrality bins are 0%-20%, 20%- 
40% and 40%-93%, and the data points are placed at the 
mean Apart calculated from a Glauber simulation. The data 
points for the Gu-going direction, 1.2 < y < 2.2, are shifted 
along the x-axis to Apart-1-3 to make the points visible, while 
the Au-going direction, —2.2 < y < —1.2, remains unshifted. 
The values are shown in Table HV] 
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FIG. 6. (color online) Invariant yield as a function of trans¬ 
verse momentum for 1.2 < \y\ < 2.2 and 0%-93% centrality. 
The Pt bins are 1 < pr < 2.5 and 2.5 < pt < 5 GeV/c, 
and the data points are placed at the mean pT of the bin. 
The Cu-going direction corresponds to the forward rapidity, 
1.2 < y < 2.2, while the Au-going direction corresponds to 
the backward rapidity, —2.2 < y < —1.2. The values are 
shown in Table 13 


nuclear matter effects in this collision system, the trans¬ 
verse momentum dependence of the nuclear modification 
factor is shown in Fig. 0 The data points are placed at 
the mean pp of the bin. Here, the nuclear modification 
is calculated over integrated centrality, but it should be 
noted that the data are dominated by central collisions. 
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TABLE IV. Invariant yield as a function of centrality for 1 < pr < 5 GeV/c and 1.2 < \y\ < 2.2. The first value represents 
the statistical and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. An additional 
±3% Type-C global systematic uncertainty also applies to the yields. The last column summarizes the forward/backward ratio 
shown in Fig. 1111 The forward/backward ratio has no Type-C systematic uncertainty. 


Centrality Bin 

(A^part) 

BR^ (Cu-going) 

(Au-going) 

forward/backward ratio 

0%-20% 

20%-40% 

40%-93% 

154.8 ±4.1 

80.4 ±3.3 

19.5 ±0.5 

(7.3 ± 7.5 ± 1.1) X 10“® 

(1.2 ±0.3 ±0.2) X IQ-* 

(1.5 ±0.6 ±0.2) X 10“® 

(3.4 ± 1.0 ±0.5) X 10“'* 

(1.2 ±0.3 ±0.2) X IQ-* 

(2.7 ± 0.7 ± 0.4) X 10“® 

0.2+° ®± < 0.1 
1.0l°;t±0.1 

0.6+°;3 ±0.1 


TABLE V. Invariant yield as a function of pr for 0%-93% centrality and 1.2 < |j/| < 2.2. The first error represents the statistical 
and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. An additional ±5.8% Type-C 
global systematic uncertainty also applies. 


p“- 

p?:"" 

(Cu-going) 

(Au-going) 

(GeV/c) 

(GeV/c) 

(GeV/c)-® 

(GeV/c)-® 

1.0 

2.5 

(2.7 ±0.8 ±0.4) X 10"® 

(5.5 ±0.9 ±0.8) X 10"® 

2.5 

5.0 

(1.8 ± 1.0 ±0.3) X 10"’’ 

(4.1 ±1.0 ±0.6) X 10"’’ 


TABLE VI. Invariant yield as a function of rapidity for 0%-93% centrality and 1 < pr < 5 GeV/c. The first error represents 
the statistical and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. An additional 
±5.8% Type-C global systematic uncertainty also applies. 



l^pax 

BR^ (Cu-going) 

BR^ (Au-going) 

1.8 

2.2 

(6.4 ±3.1 ±0.9) X 10"® 

(1.1 ±0.2 ±0.2) X lO"'’ 

1.2 

1.8 

(5.3 ±2.3 ±0.8) X 10"® 

(1.1 ±0.3 ±0.2) X 10"^ 


TABLE VII. Nuclear modification factors as a function of centrality for 1 < pr < 5 GeV/c and 1.2 < \y\ < 2.2. The first error 
represents the statistical and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. An 
additional ±15% Type-C global systematic uncertainty also applies. 


Centrality Bin 

(Acoll) 

7?CuAu (Cu-going) 

RcuAu (Au-going) 

0%-20% 

313.8 ± 28.4 

0.4 ±0.4 ±0.1 

1.7 ±0.5 ±0.3 

20%-40% 

129.3 ± 12.4 

1.4 ±0.4 ±0.3 

1.4 ±0.3 ±0.3 

40%-93% 

21.6 ± 1.0 

1.1 ±0.5 ±0.2 

1.9 ±0.5 ±0.3 


TABLE VIII. Nuclear modification factors as a function of pr for 0%-93% centrality and 1.2 < |p| < 2.2. The first error 
represents the statistical and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. An 
additional ±11% Type-C global systematic uncertainty also applies. 


p?:- (GeV/c) 

p?:“ (GeV/c) 

RcuAu (Gu-going) 

RcuAu (Au-going) 

1.0 

2.5 

1.1 ±0.4 ±0.2 

2.3 ±0.4 ±0.4 

2.5 

5.0 

0.6 ±0.4 ±0.1 

1.4 ±0.3 ±0.2 


TABLE IX. Nuclear modification factors as a function of rapidity for 0%-93% centrality and 1 < pt < 5 GeV/c. The first 
error represents the statistical and Type-A systematic uncertainties, while the second is the systematic uncertainty of Type-B. 
An additional ±11% Type-C global systematic uncertainty also applies. 



|y|max 

RcuAu (Cu-going) 

RcuAu (Au-going) 

1.8 

2.2 

1.2 ±0.6 ±0.2 

2.1 ±0.5 ±0.3 

1.2 

1.8 

0.7 ±0.3 ±0.1 

1.4 ±0.4 ±0.2 
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FIG. 7. (color online) Invariant yield as a function of rapidity 
for 1 < pt < 5 GeV/c and 0%-93% centrality. The rapidity 
bins are 1.2 < \y\ < 1.8 and 1.8 < \y\ < 2.2 and the data 
points are placed at the mean y of the bin. The Cu-going 
direction covers the region 1.2 < y < 2.2, while the Au-going 
direction covers the region —2.2 < y < —1.2. The values are 
shown in Table lYD 


FIG. 9. (color online) The nuclear modification factor f?cuAu 
as a function of transverse momentum for 1.2 < \y\ < 2.2 and 
0%-93% centrality. The pr bins are 1 < pr < 2.5 and 2.5 < 
Pt < 5 GeV/c, and the data points are placed at the mean 
Pt of the bin. The Cu-going direction corresponds to the 
forward rapidity, 1.2 < y < 2.2, while the Au-going direction 
corresponds to the backward rapidity, —2.2 < y < —1.2. The 
values are shown in Table IVTIIl 



part 


FIG. 8. (color online) The nuclear modification factor 7 ?cuAu 
as a function of the number of participating nucleons for 1.2 < 
\y\ < 2.2 and 1 < pt < 5 GeV/c. The centrality bins are 0%- 
20%, 20%-40% and 40%-93%, and the data points are placed 
at the mean Apart calculated from a Glauber simulation. The 
data points for the Cu-going direction, 1.2 < y < 2.2, are 
shifted along the x-axis to Apart-|-3 to make the points visible, 
while the data points for the Au-going direction, —2.2 < y < 
— 1.2, remain unshifted. The values are shown in Table rvm 


There is an enhancement at low pp in the Au-going di¬ 
rection. In the Cu-going direction, i?cuAu is consistent 
with unity. The enhancement in the Au-going direction 
is similar in scale to that observed in the Au-going di¬ 
rection in d-|-Au collisions [l5 | , indicating similar nuclear 


modification between the two collision systems. 



FIG. 10. (color online) The nuclear modification factor AcuAu 
as a function of rapidity for 1 < pt < 5 GeV/c and 0%- 
93% centrality. The rapidity bins are 1.2 < \y\ < 1.8 and 
1.8 < IpI < 2.2 and the data points are placed at the mean 
y of the bin. The values are shown in Table nH Also in¬ 
cluded are previous PHENIX results for ij) mesons in d-|-Au 
collisions [l3l| represented by open circles and J/i/i mesons in 
Gu-|-Au collisions represented by open triangles. Positive 
rapidity, 1.2 < y < 2.2, corresponds to the Cu-going and d- 
going directions, while negative rapidity, —2.2 < y < —1.2, is 
the Au-going direction. 

Fig. [TUI shows the nuclear modification factor RcuAu as 
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FIG. 11. (color online) The forward/backward ratio as a func¬ 
tion of the number of participating nucleons for 1 < pr < 5 
GeV/c and 1.2 < \y\ < 2.2. The values are shown in 
Table HVl The Cu-going direction covers positive rapidity, 
1.2 < y < 2.2, while the Au-going direction covers negative 
rapidity, —2.2 < y < —1.2. 


a function of y for two rapidity regions, 1.2 < \y\ < 1.8 
and 1.8 < \y\ < 2.2. The data points are placed at 
the mean y of the bin. As in Fig. |9l the nuclear mod¬ 
ification factor is inclusive of centrality. The rapidity- 
dependence of i?cuAu is similar to the trend observed 
in previous (f> meson measurements in p(d)+Au colli¬ 
sions. In particular, (j) meson production is enhanced 
in the Au-going direction. None of the Cu-going points 
show significant suppression given the statistical uncer¬ 
tainties. For comparison, the PHENIX J/ij) meson re¬ 
sults in the same Cu-I-Au dataset from Ref. [T^ are 
also shown in Fig. [TUI While the closed charm shows 
suppression at both forward and backward rapidity for 
1.2 < \y\ < 2.2, the closed strangeness is enhanced at 
backward rapidity. In Cu-|-Au collisions, the J/ip me¬ 
son yield is strongly suppressed in the Au-going direc¬ 
tion compared to the (j) meson yield at the same rapidity. 
This is similar to the differences previously observed be¬ 
tween J/ip and (/) meson nuclear modification in d-|-Au 
collisions [l^. These differences could be attributed to a 
larger J/ip break up cross section, effects in the higher- 
energy-density backward-rapidity region, or changes be¬ 
tween soft and hard production mechanisms between the 
two mesons. 

The forward and backward differences can be quanti¬ 
fied by the ratio of the yield values for the forward rapid¬ 
ity (Cu-going direction) to the backward rapidity (Au- 
going direction). Fig. [TT] shows the forward/backward 
ratio as a function of participating nucleons for 1.2 < 
lyl < 2.2 and 1 < pT < 5 GeV/c. The Type-C and 
Type-B systematic uncertainties, except for the Asrec un¬ 
certainty, cancel when taking this ratio. The remaining 
systematic uncertainties are the Type-A signal extraction 


uncertainty and the Type-B Agree uncertainty. The dif¬ 
ference in suppression between the forward and backward 
rapidity is more noticeable in the most central collisions, 
0%-20%. In this centrality bin, the probability of observ¬ 
ing the forward/backward ratio greater than or equal to 
unity was found to be p-value=1.2%, corresponding to a 
statistical significance of 2.3cr. The particle multiplicity 
for central collisions should be about 20% higher in the 
Au-going direction than in the Cu-going direction [3^ . 
however, the much smaller ratio observed may indicate 
that increased recombination effects or additional ther¬ 
mal strangeness production may also occur at higher en¬ 
ergy density. In central collisions, the forward/backward 
ratio in p production (^0.2) is smaller than that in J/ip 
production (~0.8) in Cu-fAu collisions [T^ . 


V. SUMMARY 

In summary, p meson production and its nuclear mod¬ 
ification have been measured in Cu-|-Au collisions at 
200 GeV for 1.2 < \y\ < 2.2 and 1.0 < px < 5.0 
GeV/c via the dimuon decay channel. This first mea¬ 
surement of (j) meson production and its nuclear modi¬ 
fication in a heavy-ion system at forward/backward ra¬ 
pidity at RHIC extends measurements of (j) from smaller 
systems, p+p and d-|-Au, in the forward and backward 
rapidity. The invariant yields and nuclear modification 
factors have been presented here as a function of Apart, 
PT and rapidity. 

The 4> meson yields in Cu-|-Au collisions are found to 
be generally smaller in the Gu-going direction than in 
the Au-going direction. This is most pronounced in the 
most central events, 0%-20%, and at low momentum, 
1.0-2.5 GeV/c. In central collisions (0%-20%), the for¬ 
ward/backward ratio is below unity at a confidence level 
of 99%. It has been shown that these results follow a 
trend similar to what was seen previously at PHENIX in 
d-|-Au at the same rapidity and energy |13l| as well as the 
ALICE measurement in p-|-Pb collisions at larger rapidity 
(—4.46 < y < —2.96 and 2.03 < y < 3.53) and higher en¬ 
ergy (\/aTn= 5.02 TeV) [l^. While this agreement could 
imply a role for CNM effects on (p production in Cu-|-Au 
collisions, the production of </> in heavy-ion collisions for 
these kinematics is expected to have substantial contri¬ 
butions from HNM effects as which were demonstrated 
to dominate previous measurements at midrapidity for 
both Cu-|-Cu and Au-|-Au collisions [l^. A competi¬ 
tion between CNM and HNM production mechanisms 
appears relevant for p production at forward rapidity for 
heavy-ion collisions and a comprehensive description is 
needed from soft and hard physics models. Although the 
(p meson is sensitive to both CNM and HNM effects, this 
study was statistically limited, a factor that also affects 
the precise determination of the systematic uncertainties. 
A high statistics measurement and theory calculations 
are both needed in order to make conclusions about the 
various physics processes that might be at play here, in- 
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eluding modifications of strangeness production in bulk 
matter and quark recombination. 
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